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ABSTRACT 

We report a detection of the coherent distortion of faint galaxies arising from gravi- 
tational lensing by foreground structures. This "cosmic shear" is potentially the most 
direct measure of the mass power spectrum, as it is unaffected by poorly-justified as- 
sumptions made concerning the biasing of the distribution. Our detection is based 
on an initial imaging study of 14 separated 8' x 16' fields observed in good, ho- 
mogeneous conditions with the prime focus EEV CCD camera of the 4.2m William 
Herschel Telescope. We detect an rms shear of 1.6% in 8' x 8' cells, with a signifi- 
cance of 3.4cr. We carefully justify this detection by quantifying various systematic 
effects and carrying out extensive simulations of the recovery of the shear signal from 
artificial images defined according to measured instrument characteristics. We also 
verify our detection by computing the cross-correlation between the shear in adja- 
cent cells. Including (gaussian) cosmic variance, we measure the shear variance to be 
(0.016) 2 ± (0.012) 2 ± (0.006) 2 , where these Ict errors correspond to statistical and sys- 
tematic uncertainties, respectively. Our measurements are consistent with the predic- 
tions of cluster-normalised CDM models (within la) but a COBE-normalised SCDM 
model is ruled out at the 3. Oct level. For the currently-favoured ACDM model (with 
Q m = 0.3), our measurement provides a normalisation of the mass power spectrum 
of as = 1-5 ± 0.5, fully consistent with that derived from cluster abundances. Our 
result demonstrates that ground-based telescopes can, with adequate care, be used 
to constrain the mass power spectrum on various scales. The present results are lim- 
ited mainly by cosmic variance, which can be overcome in the near future with more 
observations. 

Key words: cosmology: observations - gravitational lensing, large-scale structure of 
Universe. 



1 INTRODUCTION 

Determining the large scale distribution of matter remains 
a major goal of modern cosmology. Comparisons between 
theory and observations are hampered fundamentally by the 
fact that the former is concerned with dark matter whereas 
the latter usually probes luminous matter, particularly when 
the distribution is probed by galaxies and clusters. By con- 
trast, gravitational lensing provides direct information con- 
cerning the total mass distribution, independently of its 
state and nature. As a result, lensing has had considerable 
impact in studies of cluster mass distributions (see reviews 
by Fort & Mellier 1994, Schneider 1996) and observational 
limits have improved significantly. Weak shear has now been 
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detected >1.5 Mpc from the centre of the cluster 00024+16 
(Bonnet et al. 1994), and in a supercluster (Kaiser et al. 
1998). 

Weak lensing by large-scale structure also produces 
small coherent distortions in the images of distant field 
galaxies (see Mellier 1999; Kaiser 1999; Bartelmann & 
Schneider 1999 for recent reviews). A measurement of this 
effect on various scales (defined as 'cosmic shear') would pro- 
vide invaluable cosmological information (Kaiser 1992; Jain 
& Seljak 1997; Kamionkowski et al. 1997; Kaiser 1998; Hu 
& Tegmark 1998; Van Waerbeke et al. 1998). In particular, 
it would yield a direct measure of the power spectrum of 
density fluctuations along the line of sight and thus provide 
an independent constraint on large scale structure models 
and cosmological parameters. 

Because of its small amplitude (a few percent on arcmin 
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scales for favoured CDM models), cosmic shear has however 
been difficult to detect. In a pioneering paper, Mould et al. 
(1994) attempted to detect the coherent distortion of R ~26 
field galaxies over a 9.6 arcmin diameter field and found an 
upper limit quoted in terms of the rms shear at the 4% level. 
A search for this effect is the object of active observational 
effort (Van Waerbeke et al. 1998; Refregier et al. 1998; Seitz 
et al. 1998; Rhodes et al. 1999; Kaiser 1999). At present 
however, no unambiguous detections of cosmic shear have 
been reported (see however the limited results of Villumsen 
1995; Schneider et al. 1998). 

A fundamental limitation of narrow field imaging as a 
probe of cosmic shear is that arising from cosmic variance, 
i.e. the fluctuation in the lensing signal measured with a 
limited number of pencil beam sight lines. Only through the 
analysis of image fields in many statistically-independent 
directions can this variance be overcome. Prior to such a 
measurement, it is important to demonstrate a reliable de- 
tection strategy, particularly in the presence of significant 
instrumental and other systematic effects. 

In this paper, we report the detection of a cosmic shear 
signal with 14 separated 16' x 8' fields observed with the 4.2m 
William Herschel Telescope (WHT). We provide a detailed 
treatment of systematic effects and of the shear measure- 
ment method. We test our results with numerical simula- 
tions of lensed images and quantify both our statistical and 
systematic errors. We discuss the consequence of our mea- 
surement for the normalisation of the mass power spectrum. 
Subsequent papers will extend this technique to a larger 
number of fields, reducing the limitations caused by cosmic 
variance. 

This paper is organised as follows. In we introduce 
the theory of weak lensing in the context of a cosmic shear 
survey. In §^ we discuss our observational strategy for de- 
tecting it and describe our observations taken at the WHT 
and the routine aspects of data reduction. In §Q we describe 
the generation of the object catalogue and how the image 
parameters were measured. In §^|we discuss and characterise 
distortions introduced by the telescope optics. In §[] we dis- 
cuss the point spread function and present our shear mea- 
surement method, alongside an important comparison with 
the same analysis conducted with simulated data (§Q). In 
§|i| we describe the estimator used for measuring the shear 
variance and the cross-correlation between adjacent cells. In 
we present our results. Our conclusions are summarised 
in bEoL 



2 THEORY 

2.1 Distortion Matrix 

Gravitational lensing by large scale structure produces dis- 
tortions in the image of background galaxies (see Mellier 
1999; Kaiser 1999; Bartelmann & Schneider 1999 for recent 
reviews). These distortions are weak (about 1%) and can be 
fully characterised by the distortion matrix 
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where S8i(8) is the displacement vector produced by lensing 
on the sky. The convergence k describes overall dilations 



and contractions. The shear 71 (72) describes stretches and 
compressions along (at 45° from) the x-axis. 

The distortion matrix is directly related to the matter 
density fluctuations along the line of sight by 
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where $ is the Newtonian potential, \ is the comoving dis- 
tance, Xh is the comoving distance to the horizon, and di is 
the comoving derivative perpendicular to the line of sight. 
The radial weight function g(x) is given by 
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where r is the comoving angular diameter distance, and n(x) 
is the probability of finding a galaxy at comoving distance 
X and is normalised as J dx n (x) = 1. If the galaxies all lie 
at a single distance Xs, n (x) — $(x ~ Xs) an d 
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In practice, this approximation is accurate to within 
10%, if Xs is set to the median distance of the galaxy sample. 
This is adequate given the median redshift of our galaxy 
sample is itself uncertain by about 25% (see §3.2), yielding 
an uncertainty in the predicted rms shear of about 20% (see 
Eq. [d| below). 

2.2 Power Spectrum 

The amplitude of the cosmic shear can be quantified statisti- 
cally by computing its 2-dimensional power spectrum (Jain 
& Seljak 1997; Kamionkowski et al. 1997; Schneider et al. 
1997; Kaiser 1998) . For this purpose, we consider the Fourier 
transform of the shear field 
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The shear power spectrum is defined by 

(7i(I)7j(l')) = (2»r) 3 * W (l-l')C* j (6) 

where <5' 2 ' is the 2-dimensional Dirac-delta function, and 
the brackets denote an ensemble average. It is also useful to 
define the scalar power spectrum Ci = Cj 11 + Cf 2 for the 
shear amplitude by 

(7i(l)7i(0> = (2T)V 2 J(l-l')a, (7) 

where the summation convention was used. 

Applying Limber's equation in Fourier space (Kaiser 
1998) to Equation (^|) and using the Poisson equation, wc 
can express the shear power spectrum Ci in terms of the 
3-dimensional power spectrum P(k, x) of the mass fluctua- 
tions Sp/p and obtain 

i6 {-) n ™ L dx 
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where a is the expansion parameter, and Ho and fi m are 
the present value of the Hubble constant and matter density 
parameter, respectively. After noting that Ci is also equal to 
the power spectrum of the convergence «, we find that this 
expression agrees with that of Schneider et al. (1997). The 
component-wise power spectrum is given by 
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Figure 1. Shear power spectrum for each cosmological model and 
for sources at z s = 1. Note that the SCDM spectrum is larger due 
to its higher normalisation. 



(9) 



where Ui(X) — {cos(2A), sin(2A)} and A is the angle of the 
vector 1, counter-clockwise from the Zi-axis. 

A measurement of the power spectrum enables differ- 
entiation between the different cosmological models listed 
in Table [j] Standard Cold Dark Matter (SCDM) is ap- 
proximately COBE- normalised (Bunn & White 1997) , while 
the other variants are approximately cluster normalised 
(crs^m 3 = 0.6 ± 0.1; Viana & Liddle 1996). For each model 
we compute the non-linear power spectrum using the fitting 
formula of Peacock & Dodds (1996). The resulting power 
spectra are shown in Figure |l| for sources observed at z s = 1. 

2.3 Cell- Averaged Statistics 

For our measurement, we will consider statistics of the shear 
averaged over angular cells on the sky. This has the advan- 
tage of diminishing the impact of systematic effects (Rhodes 
et al. 1999) and allows extension in later surveys to minimise 
cosmic variance. The average shear 7 y i in a cell can be writ- 
ten as 
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where W{9) is the cell window function and is normalised 
as / d 2 9 W{9) = 1. 

It is convenient to define the Fourier transform of the 
window function as 



Wi = I d 2 9 W{9)e ae . 



For a square cell of side a, this is 
sin(aii) \ / sin(a/2) 
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To a good approximation, we can ignore the small azimuthal 
dependence of the window function and approximate 



Wi 



al/s/2 J 



(13) 



Let us consider 2 cells separated by an angle 9. We are 
interested in the correlation function 
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As is the case in our experiment, we take the separation 
vector 9 to lie along the #i-axis (or equivalently along the #2- 
axis). By taking Fourier transforms and using Equation (^), 
we thus obtain 
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As noted above, we have ignored the azimuthal dependence 
of the window function Wi. In particular, the shear variance 
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We will denote the component-wise covariances between two 
adjacent cells by 
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and their modulus by a\ = cr^i + °"x2- The values of these 
statistics for each model are listed in Table ^ for our cell 
size of a — 8'. The rms shear is of the order of 1% for the 
cluster-normalised models and of about 2% for the COBE- 
normalised model. The cross-correlation rms is about half 
the zero-lagvalue (c.f. Schneider et al. 1997). 

Figure ti shows the dependence of <r 7 on the source red- 
shift z a and as for the ACDM model (again for a = 8'). 
The range chosen approximately reflects the likely uncer- 
tainty in these parameters for our experiment. Importantly, 
the rms shear is more sensitive to as. A 10% uncertainty in 
the source redshift results in a 8% uncertainty in ct 7 . For 
this model, the dependence of <r 7 is very well approximated 
by 



0.81 1.25 



(18) 

Seljak (1997). 



<r 7 ~ 0.01152, 

in agreement with the scaling laws of Jain 
3 DATA 

3.1 Survey Strategy 

In order to detect and ultimately measure the cosmic shear, 
an array of deep imaging fields is required. These must be 
randomly placed on the sky to provide a fair sample, and 
should be well separated in order to be statistically inde- 
pendent, from the point of view of cosmic variance. As men- 
tioned in it is expedient to distinguish between a de- 
tection based on a careful analysis of a few fields, noting 
carefully the systematic effects, before embarking upon an 
exhaustive measurement survey utilising a larger number of 
fields to beat down the uncertainties arising from cosmic 
variance. With these factors in mind, we now discuss our 
strategy and observations using the William Herschel Tele- 
scope (WHT). 
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Table 1. Cell-averaged statistics for each cosmological model (with z s = 1) 



Model 


£7 m 


n A 


erg r 


<T 7 (%) 


(%) 


<7*1 (%) 


^x2 (%) 


SCDM 


1.0 





1 0.50 


2.60 


1.62 


1.23 


1.05 


tCDM 


1.0 





0.6 0.25 


1.25 


0.86 


0.64 


0.58 


ACDM 


0.3 


0.7 


1 0.25 


1.15 


0.71 


0.54 


0.46 


OCDM 


0.3 





1 0.25 


1.04 


0.62 


0.48 


0.39 



0.016 
0.014 
j- 0.012 
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Figure 2. Dependence of the rms shear on the source redshift z B 
and the power spectrum normalisation erg . The cell was chosen to 
be a square of side a = 8'. 



A bank of appropriate fields were selected for obser- 
vation with the WHT prime focus CCD Camera (field of 
view 8' x 16', pixel size 0.237", EEV CCD) in the R band. 
This photometric band offers the deepest imaging for a given 
exposure time with minimal fringing. Fields were selected 
using the Digital Sky Survey by choosing coordinates ran- 
domly within the range appropriate for the time of observa- 
tions. Each field was retrospectively checked to see whether 
it contained large galaxies (^> 5 arcsec) (which would oc- 
cult a significant fraction of the imaging field) or promi- 
nent groups/clusters (located using the NASA/IPAC Ex- 
tragalactic Database) on a scale comparable to that under 
study (~8 arcmin). There is, of course, a danger of over- 
compensating by exclusion in this respect but, fortunately, 
none of the originally-chosen fields were discarded according 
to the above criteria. 

The fields were further required to be > 5° away from 
one other, in order to ensure statistical independence (c.f. 
figure jj], where the power is small for I < 10). 

Using the APM and GCC catalogues, we ensured that 
the fields contained no stars with R < 11 (in order to avoid 
large areas of saturation and ghost images). On the other 
hand, we required the fields to contain ~ 200 stars with 
R < 22 in order to map carefully the anisotropic PSF and 
the camera distortion across the field of view. In order to 
achieve this, the fields were chosen to be at intermediate 
Galactic latitudes (30° < b < 70°; see Table |^). A calibration 



of the stellar density at limits fainter than the APM and 
GCC catalogues was obtained from a test WHT image (see 
below) . 

The final constraint on field position was our desire to 
observe each field within 20° of the telescope's zenith dur- 
ing the observing run; this reduces image distortion intro- 
duced by telescope and instrument flexure. This criterion 
was relaxed for the fields VLT1, CIRSI1 and CIRSI2 (see 
nomenclature below). 

Table |2| summarises the positions and Galactic latitude 
of the fields which are used in this paper. Two fields are 
in common with the VLT (Mellier et al., in preparation) 
and HST STIS (Seitz et al. 1998) cosmic shear programmes, 
allowing future comparisons with these programmes. A fur- 
ther two fields spanned the Groth Strip (Groth et al. 1998; 
Rhodes 1999) a deep survey conducted with HST, which has 
previously been studied for cosmic shear detection (Rhodes 
et al. 1999) . Finally, two fields were chosen to be in common 
with the current CIRSI photometric redshift survey (Firth 
et al., in preparation) to give us clearer understanding of the 
redshift distribution of objects in our fields at a later date. 

An exposure time of 1 hour on the WHT enables the de- 
tection of R=26 objects with a signal-to-noise of 5.8 in 0.8" 
seeing. This limit should correspond to a median redshift of 
about z s ~1.2. In our eventual analysis, we will introduce a 
brighter limit so as to keep only resolved galaxies (referred 
to as the survey sample). This serves to reduce the median 
redshift to about 0.8 (see §9.3). We however note from fig- 
ure ^| that our expected shear signal is not very sensitive to 
median redshift (er 7 oc z^' 8 ). In we will show that the 
resulting depth is still sufficient to detect the lensing signal. 

3.2 Observations 

We observed 14 selected fields with the WHT during the 
nights of 13-16 May 1999. For each field, a total of 4 expo- 
sures in R, each of 900s, was taken. All fields were observed 
as they passed through the meridian. 

Each exposure on a given field was offset by 10" from its 
predecessor in order to remove cosmetic defects and cosmic 
rays, and to measure the optical distortion of the telescope 
and camera (see §[| and fQ). All but two of the fields were 
observed with the long axis of the CCD pointing East- West; 
the exception being the two Groth fields, for which a 45° 
rotation (i.e. North- West orientation) was effected to align 
the WHT exposures with the HST survey (Groth et al. 1998; 
Rhodes 1999). Bias frames and sky flats were taken at the 
beginning and end of each night, and standard star obser- 
vations were interspersed with the science exposures. The 
median seeing on our used exposures is 0.81"; one exposure 
with seeing> 1.2" was excluded. 

Table H lists the seeing and 'imcat' magnitudes corre- 
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Table 2. Field Coordinates (equinox 2000) and Properties 



Field name 


RA (h:m:s) 


Dec (d:m:s) 


Galactic 


Seeing 


Magnitude 


Median 


No. Survey Galaxies 








latitude 


(arcsec) 


limit 


magnitude 


(8' X 16' field) 








(deg) 




(imcat 5<t) 


of survey 
galaxies 




WHTO 


02:03:09.31 


11:30:20.0 


-47.6 


0.59 


26.2 


23.1 


1550 


WHT3 


14:00:15.00 


10:13:40.0 


66.6 


0.82 


26.2 


23.3 


2141 


WHT5 


14:50:46.67 


20:18:03.2 


61.9 


0.76 


26.5 


23.7 


2181 


WHT7 


15:13:40.86 


36:31:30.8 


58.6 


0.83 


25.9 


23.0 


1354 


WHT11 


16:31:44.28 


27:56:30.0 


41.6 


0.85 


26.0 


23.3 


1379 


WHT12 


16:37:20.00 


20:46:30.0 


38.4 


0.90 


26.0 


23.3 


1855 


WHT1 A 

VV J- J- J- J-'rt 


1 -1 ^ ^8 

lU.i>1.10.00 




36.8 


0.99 


25 9 


23.2 


1701 


WHT16 


17:13:40.00 


38:39:19.0 


34.9 


0.78 


25.8 


23.4 


2074 


WHT17 


14:24:38.10 


22:54:01.0 


68.5 


0.63 


27.3 


24.5 


2287 


VLT1 


12:28:18.50 


02:10:05.0 


64.4 


0.71 


26.4 


23.6 


1721 


VLT2 


15:28:43.00 


10:14:20.0 


49.3 


0.79 


26.1 


23.4 


2093 


CIRSI1 


12:05:35.01 


-07:43:00.0 


60.1 


1.14 


25.4 


22.6 


1192 


CIRSI2 


15:23:37.00 


00:15:00.0 


60.4 


0.76 


26.3 


23.5 


1824 


GROTH1 


14:17:18.74 


52:20:18.5 


53.4 


0.78 


26.1 


23.4 


2237 


GROTH2 


14:15:35.00 


52:08:48.0 


44.7 


0.89 


26.1 


23.6 


1195 



sponding to 5<r detections for each field. An 'imcat' signal-to- 
noise (see Section 7) of 5.0 corresponds to a median R—26.1; 
the median magnitude of galaxies on a field is _R=25.2. To 
measure the shear, a number of cuts have to be applied to 
our object catalog (see ^). Table |^ also lists the median 
magnitude of our final sample. At our final subsample limit 
of R=23A, the median redshift is ~ 0.8 ± 0.2 (Cohen et al 
2000, §9.3). The median number density of adopted survey 
sources is 14.3 arcmin -2 (see section ^). 

In addition to the 13 useful fields, we had already ob- 
tained a test field (WHTO) in service time, and were also 
kindly given access to a suitable archival field, WHT17. Both 
were taken in good conditions: WHTO is a 1 hour exposure 
in the / band, whereas WHT17 is a 1.5 hour exposure in R 
(chosen to include a known quasar). Removing these fields 
does not significantly alter our results. In terms of unifor- 
mity, apart from the deeper WHT17 field, the standard de- 
viation in limiting magnitude is ~0.2 mag which we consider 
acceptable for our survey. The error on magnitude zero point 
derived from standard stars is at the much lower level of 0.03 
mag. 



3.3 Data Reduction 

The reduction of these deep images proceeded along a stan- 
dard route. A median-combined bias frame was subtracted 
from the skyflats and science exposures, and all such ex- 
posures were divided by a median unit-normalised sky- flat. 
Although the survey exposures were undertaken in the R 
band to avoid fringing, fringing is still detected at a 0.5% sky 
level. In order to remove these fringes, which could poten- 
tially introduce structure into the image ellipticities, all long 
dithered exposures for a given night (> 15 exposures per 
night) were stacked without offsetting with a sigma-clipping 
algorithm. This results in a fringe frame mapping the back- 
ground fringes but devoid of foreground objects. The fringe 
frame for the relevant night was then subtracted from each 
science exposure individually, subtracting off the multiple 
of the fringe frame found to minimise the rms background 



noise. After applying this technique, the fringes are entirely 
imperceptible, any residue having an amplitude within the 
sky background noise. We experimented with automated 
and hand- subtraction of the fringes and verified that this 
had no noticeable effect on our shear analysis. 

The mean linear astrometric offset (in fractional num- 
ber of pixels) between the four exposures was found by pro- 
ducing SExtractor (Bertin & Arnoults 1996) catalogues for 
each exposure, containing typically 2000-3000 objects. We 
used the mean offsets of the matched objects to align the 
fields. The images were shifted by the corresponding non- 
integer number of pixels using IRAF's imshift routine, tak- 
ing linear combinations of neighbouring pixels to effect the 
non-integer pixel shifts. As discussed in §Bl we find no need 
to rotate the exposures with respect to each other, or to 
make further astrometric distortions to compensate for the 
optical distortion of the instrument. 

The resulting four exposures for each science field were 
stacked with sigma-clipping. Since each exposure is 10" away 
from the others, bad columns and cosmic rays were rejected. 
The images were examined visually and remaining defective 
pixels (e.g. a star just outside field of view leading to light 
leakage onto an area of the CCD; or highly saturated stars) 
were flagged as potentially unreliable. 



4 IMAGE ANALYSIS 

We are now ready to measure the ellipticities of the galaxies 
on each field, and to apply the necessary corrections in order 
to take into account the smearing effect of the atmosphere 
('seeing') plus tracking and other instrumental distortions 
introduced by the telescope and camera optics. Only then 
can we ascertain the true cosmic shear by averaging the 
ellipticity distributions of the corrected galaxies. If no shear 
were present on a given field, the mean ellipticity would be 
zero, within the noise expected from the non-circularity of 
galaxies and pixelisation effects. If a shear is present, the 
mean ellipticity will be significant, especially when results 
are combined from many fields. 
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Figure 3. Example reduced image (CIRSI2); the field of view is 
8'x 16'. Note that in our analysis, we divide each such field into 
two 8' x 8' cells. 



A number of methods have recently been proposed to 
derive the shear from galaxy shapes (Kaiser et al 1995; 
Rhodes et al. 1999; Kuijken 1999; Kaiser 1999). Here we 
choose the most documented method, namely the KSB for- 
malism proposed by Kaiser et al (1995) and further devel- 
oped in Luppino & Kaiser (1997) and Hoekstra et al. (1998). 
While this method is known to have a number of shortcom- 
ings (Rhodes et al. 1999; Kuijken 1999; Kaiser 1999), it is 
nevertheless the simplest and is readily available. As we will 
show in §[?] using simulations, the method is suitable for our 
purposes, after a number of precautions are taken (see Bacon 
et al. 2000 for more details). We therefore use this method as 
provided by the imcat software, a numerical implementation 
of Kaiser et al. 1995. 

The first task in this process is to detect all objects 
present on the fields down to the background noise level, 
and to measure their shapes. We then wish to measure 
their polarisabilities i.e. measures of how each is affected by 
an isotropic smear (due principally to the atmosphere), an 
anisotropic smear (due to tracking errors at the telescope 



and local coaddition errors due to astrometric distortion) 
and shear (both the real gravitational shear and optical dis- 
tortions due to the telescope and camera optics) . One should 
note the distinction between smear and shear: a smear is a 
convolution of the image with a kernel, whereas a shear is a 
stretching of the image which conserves surface brightness. 
We will now describe the method for finding objects, and 
for measuring their ellipticities and their shear and smear 
polarisabilities. 



4.1 Object Detection 

For the purpose of detecting cosmic shear, it is expedient to 
divide each of our fields into 2 8' x 8' cells, since the signal 
is stronger on smaller scales (see Figure hj). Furthermore, 
the mean shear correlation between two adjacent cells is 
expected to be about 0.7% (see Section 2), and can thus 
be used to independently verify our results. 

We use the imcat software to find objects in each cell, 
and to measure their ellipticities, radii, magnitudes, and po- 
larisabilities. The hf indpeaks routine convolves the cell with 
Mexican hat functions of varying size, and maximally signif- 
icant peaks in surface brightness after convolution are desig- 
nated objects. The radius of the hat giving the largest signal 
to noise v for a given galaxy is attributed to that galaxy as 
its filter radius r g . The local sky background is estimated by 
the getsky routine, and aperture photometry is carried out 
on the objects, determining magnitude and half-light radius 
Th for all objects using the apphot routine. 



4.2 Shape Measurement 

Using the getshapes routine, we then measure the weighted 
quadrupole moments of each object which are defined as 



d x w(x)xiXjI(x) 



(19) 



where I is the surface brightness of the object, x is angular- 
distance from object centre, and w(x) is a Gaussian weight 
function of scale length equal to r g . In this fashion we obtain 
ellipticity components 



e, = Ii/T, 

where 

h = hi - 



2/ 2 i, T = hi+h 



(20) 



(21) 



We can further define e = (ef+el) 5 , where ei = e cos 2</> and 
e2 = esin2^>, where <j> is the position angle associated with 
the elongation direction of the object (anticlockwise from 
x-axis). The trace T of the quadrupole moments provides a 
measure for the rms radius d of the object, which we define 
as 



-(In +/ 22 )//o, 



(22) 



where Io = J d 2 x w(x)I(x) is the flux of the object. 



4.3 Polarisability 

The imcat software also enables us to calculate the smear 
and shear polarisabilities. In the following, we briefly review 
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their function. It is possible (see e.g. KSB 95 Appendix) to 
calculate the effects of anisotropic smearing, by replacing 
the image I(x) in (^) with a convolved (i.e. anisotropically 
smeared) image I'(x) and by finding the effect on the origi- 
nal ei. It is found that the galaxy ellipticity e 
corrected for the smear as 



smeared CEm be 



P 9 V 



(23) 



where the ellipticities are understood to denote the rele- 
vant 2-component spinor ei, and p is a measure of PSF 
anisotropy. The tensor P/ m is the smear polarisability, a 
2x2 matrix with components involving various moments of 
surface brightness. Since for stars e* orrectcd = 0, we can set 
V = (PL 



e s *mearcd, and find 

9 _ p9 ( p* 

'-smeared 1 am \* am 



(24) 



In this fashion, we can correct a galaxy ellipticity for the 

effect of anisotropic smearing, using the smear polarisability 
pa 

A sm • 

In a similar manner, we can calculate the effect of a 
shear, however it is induced. Replacing the image I(x) in 
( p^[ ) with a weakly sheared image, we find that 

e shearcd ~ initial ^~ ^shYl (25) 

where 7 denotes the two component shear (Eq. JS^]), and P 9 h 
is the shear polarisability, a 2 x 2 matrix with components 
involving various moments of surface brightness (different 
from P/ m above). 

In practice, the lensing shear takes effect before the 
circular smearing of the PSF. Luppino and Kaiser (1997) 
showed that the pre-smear shear 7 averaged over a field can 
be recovered using 



<^7> = <^ 
where 



Here, e 9 



i) 



Pah 
P* 



(26) 



(27) 



is the galaxy ellipticity corrected for smear, 



polarisabilities calculated for a star interpolated to the po- 
sition of the galaxy in question. The interpretation of the 
division in this equation is a matter of debate; our adopted 
procedure will be found in Section 7. With the smear and 
shear polarisabilities calculated by imcat, we can therefore 
find an estimator for the mean shear in a given cell. 

In summary, we can derive a catalogue of objects on a 
cell. For every object, we determine its centroid, magnitude, 
half-light and filter radii, ellipticity components and polaris- 
abilities as defined above. We can now use these catalogues 
to understand and correct for systematic effects, particularly 
for instrumental distortion and PSF-induced effects. 



our WHT fields and compare our measured distortion field 
to that predicted by the WHT Prime Focus manual (Carter 
& Bridges 1995) . We then show how the coadding smear can 
be computed from the astrometric frames. We finally quan- 
tify the impact of these effects on our lensing measurement. 



5.1 Measurement of the Astrometric Distortion 

The distortion field introduced by the telescope and camera 
optics can be measured from the astrometric shifts of objects 
observed in several frames offset by known amounts. Let x 
be the true position of an object. Let x-' be its position 
observed in frame /, without any correction for the camera 
distortion. The observed position can be written as 



= x + <5x(x — xf) 



(28) 



where <5x is the displacement produced by the distortion. 
The vector x^ is the position of the centre of frame /, and 
can be measured as the average position of all the objects 
found in the image. We assume that the displacement field 
<5x is the same for all frames. 

The position of this object observed in another frame 
/' is x^ = x + <5x(x — ). Here, x^ is the centre position 
of the new frame, which is assumed to be displaced from 
frame / only by a translation. (This formalism can be easily 
extended to include a rotation of the frames about their 
centre, but this effect is negligible in our case). If the offset 
x^ — x-^ is small compared to the scale on which 5x varies, 
we can expand this last expression in Taylor series and get 



*(x / - x^ 



where 



d(Sxj) 
dxj 



(29) 



(30) 



is the distortion matrix at the location of the object as de- 
fined in Equation (|l|). Following the lensing conventions, the 
distortion matrix can be parametrised as 



K + 71 72 + p 
72 — P K — 71 



(31) 



where k and 7^ are the spurious convergence and shear intro- 
duced by the geometrical distortion. We have included the 
rotation parameter p, which, unlike the case of gravitational 
lensing, does not necessarily vanish. 

The 4 free parameters of the distortion matrix can thus 
be measured from the position of an object in 3 frames /, /' 
and /". This can be done by solving the system of 4 indepen- 
dent equations formed by equation ( |29| ) and its counterpart 
for / and /". The system will not be degenerate, if the offsets 
5j/ _ a nd jj/ _ are no £ p ara u e L 



5 INSTRUMENTAL DISTORTIONS 

The instrumental distortion induced by the optical system 
of the telescope must be accounted for. If left uncorrected, 
this effect can indeed produce both a spurious shear and 
a smearing during the coadding process. In the following, 
we first present our method to measure the distortion using 
dithered astrometric frames. We then apply this method to 



5.2 Distortion Field for the WHT Prime Focus 

First we can compute the expected instrumental distortion 
using the specifications in the WHT Prime Focus manual 
(Carter & Bridges 1995). The displacement field is expected 
to be radial with an amplitude of Sx = ar 3 f , where r is the 
distance from the optical axis (located at (1076.13, 2010.7) 
pixels), f = r/r is the associated unit radial vector, and a ~ 
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4.27x 10 10 pixels 2 . Using this expression in Equations ( ^p| - 
|3l| ), we can compute the distortion parameters to be 

K = 2ar 2 

2 «r 

7, = ar e z 
P = 0, 



0.1% Shear 



(32) 



where e[ = {r\ — rf, 2r\r-z\ j (r\ + r 2 ) is the unit radial ellip- 
ticity vector. This therefore predicts a radial instrumental 
shear with an amplitude growing like r 2 , reaching 7 ~ 0.001 
at the edge of the chip. This expected shear pattern is shown 
on figure || 

Figure ^] also shows a typical instrumental shear pat- 
tern measured in one of our fields. This was derived using 
the method described above applied to 3 astrometric frames 
dithered by about 10" and containing about 15 objects per 
square arcmin. The uncertainty for the mean shear compo- 
nent ji in each of the 2' x 2' cell is of about 0.0005. Astromet- 
ric measurements thus allow us to measure the instrumental 
distortion with very high accuracy. 

The measured shear pattern is also approximately ra- 
dial and agrees well with the expected pattern. More impor- 
tantly, it also has an amplitude of at most 0.001 throughout 
the field. We have inspected all of our fields in this manner, 
and have found only small field-to-field variations (of about 
0.002) for the shear patterns. In all fields, the maximum in- 
strumental shear is only 0.003 in single 2' x 2' cells. This 
number would be even smaller, for an average over a larger 
area. We also compared the convergence k and p patterns to 
that expected from the WHT manual (Eq [p2|). We again 
found good agreement with small field-to-field variations of 
about 0.002. The origin of these variations is unknown but 
could arise perhaps from telescope flexure. For our purposes, 
however, it is quite clear that the instrumental distortion is 
much smaller than the expected lensing signal. We therefore 
neglect this component in the subsequent analysis. 



5.3 Smear Arising from Co-Addition 

If left uncorrected, instrumental distortions can also pro- 
duce a systematic effect on the shapes of galaxies, during 
the coadding process. The images of a galaxy from each (dis- 
torted) frame will be slightly offset from one another, and 
will therefore combine into a blurred coadded image. Here, 
we show that this effect is equivalent to a convolution (or 
smear) by an additional kernel. Since this effect will equally 
affect the stars in the field, it will be corrected for by the 
PSF correction described in §[| It is nevertheless important 
to estimate the amplitude of this effect, and to ensure that 
it does not dominate the dispersion of the PSF anisotropy. 

Let us consider the image of an object which appears 
on Nf frames. As before, let x-* be its centre position on 
frame / (after correcting for a translation but not for the 
distortion). Let us choose the centre of our coordinate sys- 
tem to coincide with the centre-of-light x° = y\ xf /Nf of 
the coadded image. The coadded surface brightness is then 



/=1 



(33) 



where 7(x) is the (undistorted) surface brightness of the 
object, and the factor of N7 1 was added for convenience. 
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Figure 4. Expected (top) and measured (bottom) instrumental 
shear pattern for the WHT Prime Focus. The expected pattern 
was derived from the distortion model given in the WHT Prime 
Focus manual (Carter & Bridges 1995). The observed pattern was 
measured using 3 astrometric frames in one of our fields. 



Note that the effect of the distortion on the object shape 
in individual frames was treated separately in the previous 
sections, and was thus ignored in this expression. It is easy 
to see that I' can be written as a convolution of / with the 
kernel 



(2) 



(x-x>), (34) 
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where M 2 ' is the 2-dimensional Dirac-Delta function. 

To estimate the amplitude of the effect, it is conve- 
nient and sufficient to consider the normalised unweighted 
quadrupole moments 



Jij = I d x XiXjl(x) 



d xl{x) 



(35) 



(see Eq. |l^|) of the undistorted image, and similarly for the 
moments J[j of the coadded image. The unweighted mo- 
ments of the kernel Q(x) are simply 



N, 



(36) 



Because I, /' and Q are simply related by a convolution, 
their respective quadrupole moments are related by JL = 
ij (see e.g. Rhodes et al. 1999). The rms radius d' 
) of the coadded image is thus given by 




- d"„ 



(37) 



where d and d q are the rms radius for the undistorted image 
and for the kernel respectively. For simplicity, let us consider 
an object which is intrinsically circular. The ellipticity e' t of 
the coadded image (see Eq. |2fj) is then given by (ibid) 



d? + d? q 



(38) 



where ef is the ellipticity of the kernel. 

Turning to the specific case of the WHT observations, 
let us consider the ellipticity produced by the coadding 
smear on a star observed on 4 frames with a 0.7 arcsec cir- 
cular seeing. Note that the effect will be smaller for galaxies 
which are extended, and so the following estimate should 
be considered as an upper limit. For simplicity, we conser- 
vatively assume that the seeing has a gaussian profile. We 
inspected all our fields and found that the astrometric off- 
sets between the different frames was always smaller than 

0. 3 pixels. Using Equation (^) we calculated the change 
(d! — d)/d in the radius of the star which is always less 2%, 

1. e. negligible compared to intrinsic changes in the seeing 
size. Using Equation (^) we also computed the induced el- 
lipticity e of the star and found it to be of the order of 0.01 
and always less than 0.03. This is considerably smaller than 
the rms dispersion in the PSF ellipticity that we measure in 
our fields (about 0.07, see iQ) , which must therefore be due 
to other effects (tracking errors, atmospheric effects, etc). 

Again, we can conclude that smear arising via instru- 
mental distortions during image coaddition is negligible. 



In this section, we first address the anisotropic component 
of the contribution to the point spread function (PSF) and 
then the isotropic part, thus determining a measure of the 
true (corrected) mean shear in each cell. 

Although our recipe for measuring the true shear is 
straightforward, it is the success of our simulations described 
in §Q which provides justification that our results are reliable 
at the necessary 1% level. 



6.1 Anisotropic Correction 

Our approach is to use equation (|2^) to remove the 
anisotropic component of the smearing induced in the galaxy 
images. However, we must first remove the extraneous noise 
detections in our imcat object catalogue, find appropriate 
well-defined subcatalogues of stars and galaxies, and gener- 
ate a functional model for the stellar ellipticities and polar- 
isabilities over the field of view. 

Firstly, we need to remove noisy detections. We applied 
a size limit r g > 1.0 to reject the extraneous very small ob- 
ject detections which imcat finds. We also applied a signal- 
to-noise v > 15.0 cut; see § 3.2 for justification of this ap- 
parently very conservative cut. To reduce the noise in our 
measurement, we also remove highly elliptical objects with 
e > 0.5. 

Stars were identified using the non-saturated stellar lo- 
cus on the magnitude-7"h plane (see figure^), typically with 
R ~ 19-22. The distribution of stellar ellipticity over a typ- 
ical field is shown in figure |H| for this field we find e ~ 0.07 
with only slow positional variations across the field. Al- 
though the pattern varies from field to field, it is smooth 
in all cases. The rms field-to-field stellar ellipticity is rela- 
tively large, a e * — 0.068, and must therefore be removed 
with care. 

In order to use ([m]) to correct for these elongations, we 
must estimate the positional dependence of stellar ellipticity 
and polarisability by interpolation. We adopted an iterative 
approach to this problem. We first fitted a 2-D cubic to the 
measured stellar ellipticities, plotted the residual ellipticities 



Jit 



and re-fitted after the removal of extreme 



outliers (caused by galaxy contamination, blended images 
and noise). 

Figures | and | show the stellar ellipticity residual for 
the field CIRSI2. Although the mean spurious ellipticity in- 
duced by the instrumental effects is e\ ~ —0.009, e 2 — 0.052 
over the field, the residual ellipticity after correction is only 
e{ es = (0.6 ± 1.2) x 10" 3 , e r 2 es = (2.5 ± 1.2) x 10~ 3 . Despite 
the fact that the initial stellar ellipticities on our images are 



considerable (a e 



0.068), 



is thus found to be very 



6 CORRECTION FOR THE POINT SPREAD 
FUNCTION 

In order to measure the systematic alignment of faint back- 
ground galaxies due to lensing by large-scale structure, we 
need to account for more than simply the geometric dis- 
tortions discussed in the previous section. We also need to 
correct for the effect of the varying atmospheric conditions 
throughout each exposure and imperfections in telescope 
tracking, leading to an anisotropic smearing of the images. In 
addition, the isotropic smear arising from seeing circularises 
small galaxies, thereby weakening the sought-after signal. 



cess arises because of the smoothness of the variation in the 
stellar ellipticity across each field. The small residuals will 
contribute negligibly to the mean shear. 

To check the robustness of our anisotropy correction, 
we used half of the selected stars on a field, distributed 
uniformly across the field of view, to correct the PSF 
anisotropy; we then compared the final shear measurement 
obtained for this field with that obtained after anisotropy 
correction with the other half of the stars. We found that 
the final measured shear differed by only 0.1%. 

At this stage we further discard 4 of our cells for which 
our PSF interpolation model is not satisfactory. This is 
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Figure 5. Top: Stellar ellipticity distribution for the field of fig- 
ure ^ (CIRSI2). The mean value observed is e* ~ 0.07. Bottom: 
Residual stellar ellipticity after correction. The residual mean el- 
lipticity is e rea ~ 2.6 X 10" 3 . 



Figure 6. Effectiveness of the correction for PSF anisotropy as 
applied to the field CIRSI2. Stars initially have a range of elliptic- 
ities across the field of view (dots). After polynomial fitting and 
correction (see text) the stars have small mean ellipticity close to 
zero (e rea ~ 2.6 X 10 — 3 ; crosses) 

due to Tg (and consequently P s * m ) showing strong gradients 
across the cell, or due to an insufficient number of stars in 
an area of the cell leading to poor fitting of the PSF model. 

In order to correct galaxies for anisotropic smear, we not 
only need the fitted stellar ellipticity field, but also the four 
component stellar smear and shear polarisabilities as a func- 
tion of position. Here a 2-D cubic is fit for each component of 
Psm and P* h . Galaxies are then chosen from the magnitude- 
rh diagram by removing the stellar locus and objects with 
v < 15, r g < 1, e > 0.5, as described above. From our fitted 
stellar models, we then calculate e*, P s * m and P* h at each 
galaxy position, and correct the galaxies for the anisotropic 
PSF using equation (pij). As a result, we obtain ej? orrected for 
all selected galaxies in each cell. 

6.2 Isotropic Correction and Shear Measurement 

The final correction arises because of atmospheric seeing 
which induces an isotropic smear. Clearly small objects suf- 
fer more circularisation by the isotropic component of the 
smear than larger objects. The goal now is to correct for 
this bias as well as to convert from corrected ellipticities to 
a measure of the corresponding shear, using P 7 as introduced 
in section 5, in equations (^H|) to (|2^). 

We first calculate P 7 for the galaxies. We opt to treat 
P* h and Pgrn as scalars equal to half the trace of the re- 
spective matrices. This is allowable, since the non-diagonal 
elements are small and the diagonal elements are equal 
within the measurement noise (typical P£ m ,ii,22 ~ 0.10, 

fsm, 12, 21 < 5 X 10 , -P s \oO,ll ~ 1-1) Psh, 12,21 < 0.01). 

With this simplification, we calculate P 7 according to 
equation (^?]) . P 7 is typically a noisy quantity, so we fit it as 
a function of r g . We choose to treat P 7 as a scalar, since the 
information it carries is primarily a correction for the size of 
a given galaxy, regardless of its ellipticity or orientation. We 
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Figure 7. The anticorrelation of e* and 7^ plotted for all cells, 
where i=l,2 have been superposed, for (top) a 1/ > 5 cut, and 
(bottom) a i/ > 15 cut. Note the trend for v > 5. 



thus plot P* 1 and P 22 together against r g , and fit a cubic 
to the combined points. Moreover, since P 7 is unreliable for 
objects with r g measured to be less that r* g , we remove all 
such objects from our prospective galaxy catalogue. Finally, 
we calculate a shear measure for each galaxy as (c.f. Eq. 
@) 



(39) 



where the P 1 is the fitted value for the galaxy in question. 

Because of pixel noise, a few galaxies yield extreme, un- 
physical shears 7 s . To prevent these from unnecessarily dom- 
inating the analysis, we have removed galaxies with 7 s > 2. 
We then calculate the mean 7 = (7 s ) and error in the mean 



for this distribution, giving us an esti- 



mate for the mean shear in each cell and its uncertainty. 

At this point we encountered an interesting trend. We 
found that a signal/noise cut at v > 5 (as opposed to our 
conservative v > 15) reveals a strong anti-correlation be- 
tween the mean shear 7, and the mean stellar ellipticity e*. 
This can be seen clearly in figure ^. To assess the significance 
of this effect, we use the correlation coefficient 



Ci 



:><7i> 



<r(e*)<r(7i 



(40) 



For a v > 5 cut we find Ci = -0.83, C 2 = -0.80, which, 
for 32 cells, corresponds to a > 3a effect. This is clearly 
significant, and is due to an over-correction of the PSF for 
small galaxies (in Eq. However, for a cut at v > 15, 

this reduces to G\ = —0.31, C2 = —0.38, corresponding to a 
1.7 — 2.2cr significance for the correlation, which is no longer 
significant. We will take this anticorrelation into account in 
our final results. 



7 SIMULATIONS 

In order to verify our analysis, we have conducted a detailed 
study of simulated data. The principal aim is to check that 
the shear we impose on simulated images is recovered by the 
detection method described above in the context of the ac- 
tual observations. A detailed description of our simulations 
will be found in a second paper (Bacon et al. 2000). Here 
we describe the relevant details. 

We have attempted to create a realistic simulation of 
a WHT field, with appropriate counts, magnitudes, ellip- 
ticities and diameters for stars and galaxies, including the 
effects of seeing, tracking errors, pixelisation, and an input 
shear signal. 

One approach to this problem would be to directly use 
sheared Hubble Space Telescope (HST) images degraded to 
ground-based resolution. However, a test of our signal to 
the required precision requires an area which is too large to 
be available in current HST surveys. We have thus instead 
chosen a monte-carlo approach, in which large realisations of 
artificial galaxy images are drawn to reproduce the statistics 
of existing HST surveys. Specifically, we used the resolved 
image statistics from the Groth Strip, a deep HST survey 
(Ebbels 1998, Rhodes et al 1999). This HST survey sampled 
at 0.1 arcsec effectively gives us the unsmeared (i.e. before 
convolution with ground-level seeing) ellipticities and diam- 
eters of an ensemble of galaxies. The Groth Strip is a set 
of 28 contiguous pointings in V and J; it covers an area of 
approx. 108 arcmin 2 in a 3'. 5 x 44'. region. The magnitude 
limit is / ~ 26, and the strip includes about 10,000 galax- 
ies. We use a SExtractor catalogue derived from the entire 
strip by Ebbels (1998), containing magnitude, diameter, and 
ellipticity for each object. 

We model the multi-dimensional probability distribu- 
tion of galaxy properties (ellipticity - magnitude - diameter) 
sampled by this catalogue, and draw from it a catalogue of 
galaxies statistically identical to the Groth strip distribu- 
tion. We normalise to the median number density acquired 
in our observed WHT fields, and spatially distribute the 
galaxies with a uniform probability across the field of view. 
Star counts with magnitude are modelled from the WHT 
data itself, since the Groth strip does not contain enough 
stars to create a good model. 

We then shear the galaxies in our prospective simulation 
catalogue by calculating the change in the object ellipticity 
due to lensing. Here we use the relation (Rhodes et al 1999): 



(41) 



For the purposes of this paper, we ran 3 sets of simula- 
tions: the first was a null test, with zero rms shear entered for 
20 fields; the second included a 1.5% rms shear for 30 fields; 
and the third a 5% rms shear for 20 fields. This will allow 
us to check the KSB method in the weak shear-measuring 
regime. The imposed shear is uniform over a given field; this 
simplification should not affect our results, since we are pri- 
marily interested only in the mean shear measured on the 
field. We select uniform shears for each field from a Gaus- 
sian probability distribution with standard deviation equal 
to the rms shear we wish to study. 

Stellar ellipticities (i.e. tracking errors) are similarly 
chosen as uniform over a given field, taken from a Gaussian 
probability distribution with cr=0.08 (this is conservatively 
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chosen to be slightly worse than the rms stellar ellipticity of 
the stars in our data, with a— 0.07). 

We create the catalogue using the IRAF artdata pack- 
age. This takes the star and (sheared) galaxy catalogues, 
and plots the objects at the specified positions with speci- 
fied ellipticity, magnitude, diameter and morphology (only 
exponential discs and de Vaucouleurs profiles are supported; 
we input the appropriate proportion of spirals and ellipticals 
from HST morphological counts, and model irregulars as de 
Vaucouleurs profiles). 

We use the package to recreate several WHT-specific 
details: the magnitude zero point is chosen to match the 
telescope throughput, the stars and galaxies are convolved 
with the chosen elliptical PSF (seeing chosen to be 0.8"), the 
image is appropriately pixelised (0.237" per pixel), Poisson 
and read noise (3.9 electrons) are added, the appropriate 
gain (1.45 electrons / ADU) is included, and an appropriate 
sky background (10.7ADU per sec) is imposed. The PSF 
profile chosen is the Moffat profile, J(r) = (1 + (2 1//3 - 
l)(7"/r sca i c ) 2 ) _/3 , where j3 = 2.5 and r sca io is the seeing ra- 
dius; r is the distance from the centroid, transformed so that 
the profile is elliptical. This profile has wings which fall off 
more slowly than for a gaussian profile, and provides a good 
description of our seeing-dominated PSF. An example 16' x 
8' simulated field is shown in figure [| 

Once the simulated catalogues have been realised as 
images, we run these through our shear-measurement algo- 
rithm, exactly as we did for the data (see §Qand §^|). As for 
the data, we use 8'x8' cells for shear detection and measure- 
ment. Figure |^ demonstrates some of the similarity between 
the observed and simulated fields' imcat catalogues. 
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The next check is a comparison of the input shear 
for our cells against the output shear derived by the KSB 
method; our results for the 1.5% and 5% rms simulations 
are shown in figure |l(]. The figure shows that the output 
shear is clearly linearly related to the input shear, with a 
slope close to 1. As a quantitative test, we apply a linear re- 
gression fit to both components of the shear combined. For 
the 5% rms simulations we obtain 7° ut = 0.0007 + 0.847j n , 
with standard errors on the coefficients of 0.001 and 0.04, 
respectively. For the 1.5% simulations we similarly obtain 
7° ut = 0.0001 + 0.797" 1 with respective standard errors of 
0.001 and 0.091. We see that the imcat measure of shear is 
symmetrical about zero, but appears to measure the shear as 
somewhat too small; we therefore adjust our shear measures 
by dividing by 0.84 ± 0.04, and account for the uncertainty 
in our results analysis. 

For low signal-to-noise galaxies, the simulations also dis- 
play the anticorrelation betw een the ellipticities of the galax- 
ies and of the stars (see §6.2). For a v > 5 cut, the amplitude 



of this anticorrelation is consistent with that found in the 
data. As in the data, the anticorrelation is no longer signif- 
icant for a v > 15 cut. This confirms both the use of the 
simulations to test for systematic effects, and the validity of 
our signal-to-noise cut. 



Figure 8. An example simulated image (see text for details). 



8 ESTIMATOR FOR THE COSMIC SHEAR 
8.1 Shear Variance 

The amplitude of the cosmic shear can be measured by con- 
sidering the shear variance in excess of noise and systematic 
effects. In our experiment, we consider Nf 8' x 16' fields 
subdivided into iV c = 2N f 8' x 8' cells (see Figure §). Let 
7/° be the shear measured in cell c of field /. Here, c = t or 
b for the top or bottom cells in each field, respectively. This 
shear is a sum of the contributions from lensing, noise and 
residual systematic effects, and can thus be written as 



7 



fc 



lens, fc 1 noise, fc . sys, fc 

7, 1 + 7i + 7i 



(42) 



j7 lcns '- fc | 2 ) in excess of the noise variance o" 2 oisc j c ~ 
(| 7 noisc ' /c | 2 ) and systematic variance a^ fc = (|7 sy "' /c | 2 ). 
An estimator for the lensing variance is given by 



— cr„ 



lens — ^ tot ^ noise w sysi 

where the observed total variance is 



(43) 
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Figure 9. The distribution of image magnitudes and half-light 
radii, r^, for the data (top) and the simulation (bottom). These 
distributions are used for star/galaxy separation. 
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Figure 10. 7" 1 compared with 7° ut for simulated data sheared 
(top) by 1.5% rms shear; (bottom) by 5% rms shear. The dashed 
line shows the Y n = 7° ut relation; the solid line shows the best 
fit, 7 j n = 0.847°"*. 



2 

O"tot 



„/ c l 2 



(44) 



and the mean noise and systematic variances are denned by 
1 . 

No 



N c ^ 

fc 



\ ~> 2 

^noise,/c 5 



-y 



2 

"sys./c- 



(45) 



fc 



It is easy to check that this estimator is unbiased, i.e. that 



2 

*^lcns 1 



(46) 



fens if we assu me 



where the brackets denote an ensemble average 

We can also compute the variance of a 
that the variables follow a gaussian distribution. This is a 
good approximation for ^ nolsc '^ c i s i nC e we are considering 
an average over many galaxies (about 2000) in a cell. The 
systematic contribution to the shear is dominated by the 
residual anticorrelation discussed in 8 3.2 and thus has a dis- 



tribution which is close to that of the stellar ellipticities. 
The stellar ellipticities are relatively well approximated by 
a gaussian distribution. In our case, it is therefore accept- 
able to take 7i ys '^ c to be gaussian. The lensing shear ■y 1 i cns 'f c 
is however known to be non-gaussian, especially on scales 
smaller than 10' below which nonlinear density perturba- 
tions are dominant (e.g.. Jain & Seljak 1997; Gatzanaga 
& Bernardeau 1998). In principle, higher order correlation 
functions are required. These are however difficult to com- 
pute analytically on such small scales (Scoccimaro et al. 
1999; Hui 1999), and are at the limit of the resolution of 
current N-body simulations (Jain et al. 1998; Barber et al. 
1999; White et al. 1999). 

We can now compute the variance of the estimator and 

find 



1 



2 | 2 

ens ' ^noise 



2 (""xlcnsl + °"xl 



+ 



(47) 



where o" x i ons ; are the cross-correlation variances between top 
and bottom cells due to lensing (see Eq. [0). In deriving 
this expression, we have assumed that the noise and system- 
atic effects are uncorrelated from the top to the bottom cell. 
We have also used the following approximation 



1 \ ~> 2 ^ _j_ \ ~< 

f 

_L \ " 4 

N OVoise,/ 



2 

noise,/b 



(48) 



This is valid given the cells were observed in very similar 



conditions, and thus the spread in the a n 



oisc,/c 



is small. 



Terms with a 'lens' subscript in equation ( fl7| ) corre- 
spond to cosmic variance, while the other two terms corre- 
spond to the uncertainty produced by noise and systematic 
effects. If we are initially interested in a detection of cosmic 
shear, it suffices to test only the null hypothesis correspond- 
ing to the absence of lensing, i.e. to (Ji ons = Uxlensi = 0. The 
estimator variance relevant for a detection is 



> r t If 



2 

noise 



(detection). 



(49) 
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8.2 Shear Cross-Correlation 

An important aspect of our experiment is our ability to test 
the the cross-correlation between the shear measured in 2 
adjacent 8' x 8' cells (see Figure ^). As before, let 7/* and 
r y{ b be the average shear in the top and bottom portion of 
the 8' x 16' field /, respectively. The shear cross-correlation 
variance (see Eq. jlTj]) is defined by 



2 

O'xl 



(50) 



where the summation convention is used. As before, we as- 
sume that the noise and systematic effects are uncorrelated 
across the two cells. An estimator for this quantity is then 



X lens 



/ 



(51) 



It is again easy to check that it is unbiased, i.e. that 

(^Xlens) = ^Xlons (52) 



Assuming as before that the fields are gaussian, we can com- 
pute the variance of this estimator and find 



1 



2 | 2 



2 (°"xlcnsl + a xl 



+ O a 



0]. 



\2 
,) + 



(53) 



which equals 0" 2 [<7 2 ens 
null hypothesis (<7 2 lcnsl 
estimator variance for this purpose is then 



For a detection we must rule out the 
xicn S 2 = ofens = 0). The relevant 



2N f 



9 RESULTS 



+ o a 



(detection). 



(54) 



can first use this simulation to derive an independent con- 



straint on (Ts 



2 _ 2 

Osys — O"total 



— CTle 



(56) 



where we let cri cns = 0.015 i.e. the input rms shear. 

For this simulation, we find <r no i sc ~ 0.0130, Ototai — 
0.0193 (see figure ^|). The error for <j sys this time is com- 
puted as follows 



f~2~~l ^ 1 T_2 ,2,2]- 
[CTsysJ — jy [^noisc "t" ^sys ~r CJ"i cns J 



(57) 



However, since <7 t 2 otal 



°"noise < 0.015 2 , we can only find an 
upper limit for (j sys here; we find that o"[o" 2 ys ] = ±(0.0082) 2 , 
consistent with the null simulation result. Accordingly, in 
what follows, we will use the null simulation estimate for 

2 

V sys • 

Turning this around, we can use equation (^) to es- 
timate the rms shear in these simulations (ignoring our 
knowledge of the input rms shear). We obtain (using the 



null simulation estimate of <r s 



CTle 



0.013. The uncer- 



tainty in af ens is calculated using equation ( |49[ ) for detection 
and equation (^t]) for measurement. We obtain 0"[e>f CIls ] = 
(0.0060) 2 for detection, and a[af cns ] = (0.0082) 2 for mea- 
surement. Notice that this is the same value as o-[<r 2 ys ], since 
we can't independently find <r sys and o"i cns for the simula- 
tions. The measured rms shear is thus fully consistent with 
the input rms shear of 1.5% (Figure |l|). 

An analogous analysis is done for the 5% rms shear 
simulations; see table H. Again, we recover the input rms 
shear within the uncertainties. Note again that, since o" 2 ota i — 



Onoise < 0.05^, only an upper limit can be found for a ByB here. 
We can conclude that the simulations clearly show that in 
the relevant regimes, our method is unbiased. 



We now present and interpret our results, first using the sim- 
ulations, and then examining the WHT data. The following 
description is summarised in Table || for convenience. 



9.1 Simulated Fields 

We begin with the null simulations, which consists of 20 
8'x8' disjoint cells. The distribution of the shear for each 
simulated cell is shown on figure [H]. 

For the null simulation, the rms noise (Eq. Jt5|) is 
(j no i S e — 0.0103, while the observed total rms is Ototai — 
0.0113 (Eq. jlij) . The noise and total rms are indicated as 
a solid and dashed line in figure [Ti], respectively. Clearly, in 



the absence of a lensing signal, a Bys = a t , 
gives a sys ~ 0.0047. 



to that of equation^ 



we find that 



2 r 9 ^ 1 r 2 2 1 2 

a [(TsysJ ~ — [a noisc + <7 sys J 



giving G [fTg ys J 



(0.0053) 2 , so that 



In a fashion similar 



(55) 

(0.0047) 2 ± 



(0.0053) . Note that this is consistent with zero; i.e. even 
if we supposed that there were no systematics, the excess 
shear signal that we would attribute to real lensing would 
be consistent with zero. 

We can check this result against our next simulation, 
which now includes a 1.5% rms shear in 30 8'x8' cells. We 



9.2 Observed fields 

We now consider the observed fields. The distribution of 
shear for each of the 26 cells is plotted on figure [l^, along 
with circles corresponding to a no i S c and Otot. The mean 
shear components are 7 y 1 = —0.00097 ± 0.0034 , j 2 ~ 
0.0021 ±0.0034 , fully consistent with zero as they should be 
in the absence of systematic effects. What is more, we are 
measuring a total shear variance in excess of the noise. We 
now determine whether this detection is significant. 

The value for the rms noise is cr no i sc = 0.018, somewhat 
larger than in our simulations. This is due to increased noise 
from stellar ellipticity fitting and, in poorer seeing cases, 
lower number density. The total rms shear is cr t ot = 0.024, 
and using <7 sys = 0.0047 from the null simulations, we obtain 
cr lcns = 0.0156 (from Eq. ||). 

Using equation ([l9|), we find the error in <Ti ens to be 
cr[(j 2 ens ] ~ (0.0082) 2 for the statistical error only. If we also 
include the uncertainty on the systematic (by adding it in 
quadrature), we obtain c[cr? ] ~ (0.0084) 2 . We therefore 
quote our result as 



2 

^lcns 



— f lens ,mcasurcd ^ ( ^"[ ( ^"^ cns Jstatistical G\G S 

= (0.0156) 2 ± (0.0082) 2 ± (0.0047) 2 . 



(58) 



The significance of our detection of the cosmic shear is there- 
fore 
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Table 3. Shear measurement results for the simulated and WHT fields 





Sim. 


Sim. 


Sim. 


Data 




Null 


1.5% 


5% 




iV c 


20 


30 


20 


26 


"tot 


(0.0113) 2 


(0.0193) 2 


(0.0494) 2 


(0.0243) 2 


a 2 . 

noise 


(0.0103) 2 


(0.0130) 2 


(0.0102) 2 


(0.0179) 2 


2 

^sys 


(0.0047) 2 a 


c 


d 


(0.0047) 2 e 


4 

l c n & 


(0.0047) 2 b 


(0.0132) 2 


(0.0480) 2 


(0.0156) 2 


<7 2 , 

X lens 








(0.0156) 2 




(0.0053) 2 a 


(0.0082) 2 


(0.0234) 2 


(0.0053) 2 


°t°leij ( detect ) 


(0.0053) 2 


(0.0060) 2 


(0.0055) 2 


(0.0082) 2 


cr\a? 1 (measure) 


(0.0053) 2 6 


(0.0082) 2 


(0.0234) 2 


(0.0119) 2 



a assumes <7i ena =0. 
6 assumes <r ays =0. 
c assumes oi ons =0.015 
d assumes cri ens =0.05 

e Uses the null simulation value, since we can not obtain an independent estimate 



(S/AT) de tect = Zt ns =a 3.4 (59) 

In terms of measu ring the amplitude of the cosmic 
shear, we use equation mil) and find <T[0f ens ] = (0.0119) 2 ; 
including the uncertainty on the systematic we obtain 
fkiens] = (0.0121) 2 . We therefore find 

°~lcns °~lcns, measured ^ ^"[^lens] statistical i °~[°~sys] 

= (0.0156) 2 ± (0.0119) 2 ± (0.0057) 2 , (60) 

where we have included in <r sys the uncertainty in our 
KSB shear calibration (see section Q). Thus we find that 

(S/N) measure — J- . ( . 

The final measurement we can make is the cross- 
correlation covariance using equation (|5l"|). We find that 
a x i 0.0115, (7 X 2 ~ 0.0105, leading to cr x = 0.0156. 
For a detection, er[er 2 ] ~ (0.0088) 2 (Eq. |H^), so that the 
significance of the detection is (5'/A r )detect — 3.2 for the 
cross-correlation. For a measurement, o-[a- 2 ] ~ (0.0119) 2 
(Eq. H|), so (S/N) 

measure — J- . ( . 

9.3 Cosmological Implications 

A key question we must now address is the redshift distri- 
bution of our background sources. At the median magni- 
tude of the original catalogue (i?=25.2±0.2, Table 2, §3.2), 
photometric redshift estimators in various HST and ground- 
based datasets suggest a mean redshift z ~1-1.2 (Fernandez- 
Soto et al 1999, Poli et al 1999, Rhodes 1999). More impor- 
tantly, for the survey sample, which is effectively limited 
at a median magnitude of R ~ 23.4±0.2, we make use of 
the recently-completed deep spectroscopic survey of Cohen 
et al (2000) which indicates a median redshift at this limit 
z=0.8±0.2; the uncertainty here includes the observed field- 
to-field variations in this limiting depth as in Table 2. 

We can now compare these results with those predicted 
for the various cosmological models listed in Table m First, 
we compare our value of af ens = (0.016) 2 ± (0.012P (with 
errors which includes cosmic variance). We find that our re- 
sult is consistent with the cluster normalised models rCDM, 
ACDM and OCDM at the 0.6-0.9 a level, but that it is in- 
consistent with the COBE normalised SCDM model at the 
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Figure 11. Mean 71 and 72 for 20 simulated null cells. The 
dashed circle shows the noise rms, the solid circle shows the total 
rms. The difference is consistent with zero signal. 



3.0cr level. This confirms the fact that the SCDM model has 
too much power on small scales when normalised to COBE. 

The cross-correlation variance o x j cns = (0.0156) 2 ± 
(0.0119) 2 (a gain with cosmic variance included in the un- 
certainty) does not provide as strong a constraint. It is con- 
sistent with the models, with deviations of between O.lcr to 
1.4a. This results from the fact that 0" x i ens is expected to 
have a smaller amplitude than of cns in all models. It is nev- 
ertheless comforting that, within the context of the models 
considered, it is consistent with our measurement of (7 2 ens . 

We can use our measurement of of cns to constrain as, 
the normalisation of the matter power spectrum on 8 hT 1 
Mpc scales. For the ACDM model with il m — 0.3, we find 
from equation ( |lg| ) that these two quantities are related by 

a 8 = 0.894, s -°- 648 (^)°' 8 (61) 
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Figure 12. Mean 71 and 72 for 30 simulated cells with rms 1.5% 
shear. The dashed circle shows the noise rms, the solid circle shows 
the total rms. 
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Figure 13. Mean 71 and 72 for the observed cells. The dashed 
circle shows the noise rms, the solid circle shows the total rms. 



For our value of o"i ens and setting z 8 = 0.8 ± 0.2 (see §3.2) 
and propagating errors, this yields 



a 8 = 1.47 ± 0.24 ± 0.46 = 1.47 ± 0.51, 



(62) 



where the first error arises from the uncertainty in z 3 and the 
second from that of af ena . This corresponds to a 2.9a mea- 
surement of as- We can compare this with the cluster abun- 
dance determination which yield as = (0.6 ± 0.1)f^ 0,53 = 
(1.13±0.19) (773-) ' ■ We see that our result is consistent 
with this independent determination. Note that the uncer- 
tainty in z 3 does not dominate our uncertainty for a$. 



10 CONCLUSIONS 

Using 14 8' x 16' fields observed homogeneously with the 
WHT, we have detected a shear signal arising from weak 
lensing by large scale structure. Neglecting cosmic vari- 
ance (to test the null hypothesis corresponding to the ab- 
sence of lensing), we find a shear variance in 8' x 8' cells 
of (0.016) 2 ± (0.008) 2 ± (0.005) 2 , where the errors corre- 
spond to la statistical and systematic uncertainties, respec- 
tively. This corresponds to a detection significant at the 
3.4cr level. Including (gaussian) cosmic variance, the shear 
variance is (0.016) 2 ± (0.012) 2 . This is consistent with the 
value expected for cluster-normalised CDM models (ai cns = 
(1.0-1.3) x 10~ 2 ). On the other hand, the COBE-normalised 
SCDM model is rejected at the (3. Oct) level. We have verified 
our results by measuring the cross-correlation of the shear 
in adjacent cells. We find that the resulting cross-correlation 
variance for detection is (0.016) 2 ±(0.009) 2 , and for measure- 
ment is (0.016) 2 ± (0.012) 2 , in agreement with that expected 
in cluster normalised CDM models. This is consistent with 
all the models considered at the la level. 

Our measurement was derived after a careful account- 
ing of the systematic effects which can produce a spurious 
shear signal. We find that the most serious systematic effect 
is the PSF overcorrection for faint objects in the shear mea- 
surement method. We have shown however that, by keeping 
only sufficiently bright objects (S/N > 15), this effect can 
be made to be smaller than the statistical uncertainty. Our 
methods have been tested using detailed numerical simula- 
tions of the shear signal from appropriately-constructed syn- 
thetic sheared images. We find very good statistical agree- 
ment between the simulated and the observed data. An ex- 
tensive description of the simulations will be described in 
Bacon et al. (2000). 

For a given cosmological model, our measurement can 
be turned into a measurement of as, the normalisation 
of the mass power spectrum on 8 ft -1 Mpc scales. For a 
ACDM model with fl m = 0.3, we get a 8 = 1.5 ± 0.2 ± 0.5, 
where the errors are la uncertainties resulting from the 
uncertainty in the redshift of the background galaxies and 
from our measurement error, respectively. This result is con- 
sistent with the as value derived from cluster abundance 
{as = (1.13 ± 0.19) (%f)~°' 53 , Viana & Liddle 1996). 

The uncertainty in our measurement is clearly domi- 
nated by cosmic variance and statistical errors. This can be 
improved by increasing the number of fields Nf. Since the 
signal-to-noise ratio scales as \/Nf, a four fold improvement 
in Nf should yield a 6.8<r detection and a 3.4cr measurement 
of the rms shear. This, and the presence of other wide field 
cameras, offers good prospects for the improvement of the 
measurement of <78 from cosmic shear. On the other hand 
the determination of as from cluster abundance is currently 
measured only at the 6a level and is fundamentally lim- 
ited by the finite number of nearby clusters, for which ac- 
curate temperatures can be determined. In addition, it de- 
pends sensitively on the assumption of gaussian initial con- 
ditions. It is therefore likely that cosmic shear measurements 
will supplant cluster abundance for the normalisation of the 
power spectrum. With an even larger number of fields, one 
can also measure the shape of the power spectrum by looking 
at the correlation of the shear between and within nearby 
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fields. The advent of wide-field cameras will make this pos- 
sible in the near future. 
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